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ABSTRACT 

Numerical dynamo models always employ param.eter values that differ by orders of magnitude from 
the values expected in natural objects. However, such models have been successful in qualitatively 
reproducing properties of planetary and stellar dynamos. This qualitative agreement fuels the idea 
that both numerical models and astrophysical objects may operate in the same asymptotic regime 
of dynamics. This can be tested by exploring the scaling behavior of the models. For convection- 
driven incompressible spherical shell dynamos with constant material properties, scaling laws had been 
established previously that relate flow velocity and magnetic field strength to the available power. Here 
we analyze 272 direct numerical simulations using the anelastic approximation, involving also cases 
with radius-dependent magnetic, thermal and viscous diffusivities. These better represent conditions 
in gas giant planets and low-mass stars compared to Boussinesq models. Our study provides strong 
support for the hypothesis that both mean velocity and mean magnetic field strength scale as a 
function of power generated by buoyancy forces in the same way for a wide range of conditions. 
Keywords: stars: magnetic fields — stars: interiors, convection — stars: low-mass, brown dwarfs — 
methods: numerical 



1. INTRODUCTION 

Dynamo simulations aim to capture the magnetic 
field generation process in planetary and stellar interiors 
and have been very successful in quali tatively reproduc- 
ing many of the observed properti es (jBrun et al.l 120041 : 
iWicht fc TiignCTllMot [Jon^sl [20T1 . However, all such 
numerical simulations use inconsistent control parame- 
ters, either too large or too small, due to limited com- 
putational resources. As a consistency check between 
our theoretical understanding of the dynamo mechanism 
with astrophysical observations, it is of paramount im- 
portance to establish generic scaling laws which are valid 
in the relevant dynamical regimes. 

Many attempts have been made earlier to derive a 
scaling theory for the mean magnetic field, for exam 



nesq dynamos with different physical setup s and 
different bou ndary con ditions (iTakahashi et al.l 120081: 
Aubert et all J2p09: C hristensenI 120101: ISchrinner et all 



pie, on force balance considerations ([St evenson 

[199; 



1979 : 



Curtis fc Ne"si 1198^ iMizutani et all [l992; .Sand 11993 : 
Starchenko k. JonesI 120021 ). However, none of the sug- 



gested scaling laws was generally accepted. Progress 
was made in this regard due to t he increase in computa- 
tional power in the last decade. iChristensen fc AubertI 
(12006;) analyzed the results of parameter studies of dy- 
namo simulations under the Boussinesq approximation 
in a rotating spherical shell and found that mean ve- 
locity and mean magnetic field scale as a function of 
the available convective power generated via buoyancy 
forces. The power-based scaling laws (hereafter referred 
to as PBS) do not explicitly depend on rotation rate in 
case of magnetic field scaling. In the parameter range of 
the simulations, a secondary influence of diffusivities in 
the form of magnetic Prandtl number (ratio of viscosity 
and magnetic diffusivity) was found. 
Subsequent studies generalized the PBS to Boussi- 
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20121: lYadaveT all I2013D . lOlsonfc ChristensenI (12001 ) 
derived PBS for the magnetic dipole moment from 
numerical simulations and found an order of magni- 
tude agreeme nt with dynamos i n sola r sys tem planets. 
Furth ermore, IChristensen et al.) (|2009( ) and IChristensenI 
(|2010f ) observed a good agreement between PBS of mag- 
netic field from numerical simulations and observation- 
ally constrained magnetic field of Earth, Jupiter, rapidly 
rotating l ow-mass stars, and p ossibly Uranus and Nep- 
tune. Stcl zer fc Jacksonl (|2013[ ) maintain that including 
a dependence on the magnetic Prandtl number is manda- 
tory for an adequate fit, at least in the parameter ra nge 
of current numerical simulations. iDavidsonI (|2013f ) re- 
cently provided interesting theoretical arguments sup- 
porting the power based scaling laws. 

The agreement between the prediction of scaling laws 
derived from Boussinesq numerical simulations and ob- 
served magnetic fields in low-mass stars is rather puz- 
zling. Unlike the dynamo mechanism in solar type stars, 
where the strong differential rotation and magnetic field 
generation by s hear at the tachoch ne are thought to be a 
key ingredient ()Ossendriiveiil2003l) , rapidly rotating low- 
mass stars (mass < 0.35 Mq) and giant planets possibly 
harbor dynamos similar to the geodynamo, where helical 
convection columns aligned with the rotation axis are in- 
strumental. However, the hydrogen rich interiors of low- 
mass stars are vastly different from the liquid metal inte- 
riors of the Earth-like objects. The density and transport 
properties in the liquid core of the latter vary by some 
tens of percent (e .g. » 20% density chang e across Earth's 
hquid core; see iBraginskv fc Robert J |i995)) and can 
be considered constant. On the other hand, the interiors 
of gas planets and stars have significant density stratifi- 
cation and transport properties (such as electrical con- 



ductivity a nd thermal diffusivi ty) may vary by orders of 
magnitude (jFrench et al.ll201^ . 

The power based scahng laws discussed above were 
derived from a large number of Boussinesq dynamo 
simulations. Although density stratified models with 
radially-varying transport properties have been com- 
mon ly employed in the ste ll ar dy na mo commun i ty (se e 
e^gjGilman fcGlatzmaie^ (flQSlI) : iBrun et alj (|20Cil : 
[Brown ing (200^), a systematic scaling study of impor- 
tant diagnostic quantities has never been carried out. 
One of the reasons is the substantial increase in compu- 
tational requirements associated with anelastic density- 
stratified dynamo simula tions which make s parameter 
studies rather expensive (iJones et al.l I2009D . For non - 
magnetic rotating convection. i Gastine fc Wich^ (|2012j ) 
found in systematic model studies that the velocity of 
convection and of zonal flow scales in the same way for 
Boussinesq and anelastic cases. Here we extend this to 
dynamo models with density stratification. We also in- 
clude cases with different forms of variation in transport 
properties, different radial gravity profiles and different 
mechanical and magnetic boundary conditions. Our aim 
is not to model any particular class of astrophysical ob- 
ject as realistically as possible, but rather concentrate on 
generic scaling properties. Our analysis of more than 270 
numerical dynamo models shows that the same power 
based scaling laws apply to a wide variety of dynamos. 

2. EQUATIONS AND NUMERICAL SETUPS 
2.1. Anelastic MHD equations 

We consider dynamo action in spherical shells, with 
inner radius r^ and outer radius Tq, filled with an elec- 
trically conducting fluid. The aspect ratio rj is defined 
as r.i/ro- The shell rotates along a vertical axis z with 
constant angular velocity Jl. Convection inside the shell 
is driven by a fixed entropy contrast As between the in- 
ner and the outer boundary. We work in dimensionless 
units using shell thickness D — ro — ri and inverse rota- 
tion frequency Q.~^ as the fundamental length and time 
units, respectively. The density p and entropy s are non- 
dimensionalized using p{ro) = p^ and As, respectively. 
Magnetic field is scaled by ilDy/fipo, where p is the mag- 
netic permeability. 

To model the low-Mach number flows in the density 
stratifled interiors of giant planets and low-mass stars, 
we employ the anelastic approximation. It allows radial 
variation of density while filtering out the fast acous- 
tic waves out of th e system (jBraginskv &: RobertslllQQSl : 
iLantz fc FanI I1999D . The anelastic approximation as- 
sumes an adiabatic reference state, i.e. dT/dr = —g/cp, 
where g is gravity and Cp is the specific heat at constant 
pressure. The radius-dependent reference state quanti- 
ties are highlighted with a tilde on top. For the sake of 
generality we define the gravity profile as 



r r 

ffW =31 ^52-|, 



(1) 



and by using g\ or 32 appropriately we can either choose 
a linear gravity profile (51 = 1, 52 = 0), approximately 
representing a self-gravitating body with weak density 
variation, or an r~^ gravity profile [g\ = 0, (72 = 1), 
exemplifying objects with massive core. Assuming an 
ideal gas equation of state leads to a polytropic reference 



state defined by /5 = T™, where m is the polytropic index. 
Solving for T using Eq. ([!]) leads to 



T=\-Co 



^f;.-l 
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(2) 



with 



7]\t 



^(77-773) + 92(1 - ?7)' 

where Np — \i\{p{r i) j p(r o)) represents the number of 
density scale heights across the shell. 

The thermodynamic variables, density, pressure, and 
temperature are then decomposed into the sum of ref- 
erence state values and small perturbations as p + p, 
P -\- jp , and T + T respectively (iGilman fc Glatzmaien 
[19811 IBraginskv fc Roberts! [T995t ILantz fc FanI I1999D . 
The evolution of velocity u is governed by 
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where B is magnetic field and f is the radial unit vector. 
The traceless rate-of-strain tensor S is defined by 



with 



Sij = 2p 



^^j - -(5y V • U 



duj 



dx. 



(5) 
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where 5ij is the identity matrix. The entropy s evolves 
according to 



pTE 



ds 



dt 
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where the viscous heating contribution is 



Qu^2p 
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The magnetic induction is governed by 

5B 
'dt 



= V X (u X B) - ^-V X (aV X B 



p 
V • B = 0. 



(7) 



(8) 



(9) 
(10) 



In the above equations, kinematic viscosity v, thermal 
diffusivity k, and magnetic diffusivity A are normalized 
by their value at the inner boundary (r = r^); they are 

denoted by i>, k, and A. The various dimensionless con- 



trol parameters appearing in Eqs. ^l^ are: 



Rayleigh number Ra 



Ekman number E — 



Cp Vi l^i 
Vi 



Prandtl number P, 



Vi 



r — 1 

Ki 



magnetic Prandtl number P„i ~ — , 

where the "«" subscript represents values at the inner 
boundary. 

2.2. Variable properties 

Recent ah initio calculations suggest that electrical 
conductivity decreases with radius by several orders 
of m agnitude in the outer regions of Jupiter-like ob- 
jects (jFrench et al. 2012). Low-mass stars, brown dwarfs 
and massive extrasolar planets will probably show a simi- 
lar variation in electrical conductivity. An electrical con- 
ductivity profile approximately constant in the deep in- 
terior and exponentially decaying in outer portio ns rep- 
resen ts a good model for giant planets (Frenc h~et al.l 
|2012[) . We model electrical conductivity normalized by 
its value at the inner bo undary as (|G6mez-Perez et ahl 
[Mot iDuarte et al.l[20l3l) 



a(r) = 



'i+('^'»-i)(;^:-^)" ^ 



r > r„i 



(11) 



The electrical conductivity follows a polynomial in ri < 
r < r„i and decreases to a{rm) = ^m (usually 0.5) near 
rrn- The exponential decay starts for r > rm- The con- 
stant a defines the rate of decay; a is mostly equal to 
9, except of two cases where it is 1 and 25. The rel- 
ative transition radius Xm = fm/ri is 0.7, 0.8, 0.9, or 
0.95. The magnetic diffusivity A = (/i(7(r))^^ accord- 
ingly rises along the radius. We use such profiles in many 
of our anelastic dynamo simulations. 

We also simulate cases with diffusivities, viz. viscos- 
ity, thermal diffusivity, magnetic diffusivity, varying as a 
function of density: 



(zy, K, A) = {v,n,\)i. 



p{r)' 



(12) 



Such density dependence is simila r to that used i n 
many stellar dyn amo simulations (e.g. lBrun et al.l(|2004D : 
iBrowninsJ (|2008l )). 

2.3. Boundary conditions 

The mechanical boundary condition is either stress-free 
on both boundaries or rigid on the inner and stress-free 
on the outer boundary. On both boundaries constant 
entropy is imposed. We do not simulate cases with flux 
boundary conditions. For Boussinesq cases, no differ- 
ence in scaling laws had been found between fixed tem- 
perature and fixe d flux conditions (jAubert et al.l 120091 : 
lChristensenll2010l ). The inner core is either conducting 
or insulating. The magnetic field matches a diffusive so- 
lution at the inner boundary in case of a conducting inner 



core and a potential field in case of an insulating inner 
core, while it always matches a potential field on the 
outer boundary. 

2.4. Numerical technique 



The anelastic equati ons (Eqs. [BjfTOl) are time ad- 
vance d using Magic (jWichtl l200l iGastine fc Wichtl 
l2012f) . The anelastic version of MagIC has been bench- 
marked with recent comm unity-based benchmark simu- 
lations pones et al.ll20lil ). The mass-flux and the mag- 
netic field are decomposed into poloidal an toroidal parts 
as 

/5u = V X (V X Wf) + V X Xf, 

B = V X (V X yf) + V X Zf 

where W and X {Y and Z) are poloidal and toroidal 
scalar potentials for mass-flux (magnetic field). The 
scalar potentials, along with p and s, are then expanded 
using spherical harmonic functions in the azimuthal and 
the latitudinal direction. The radial expansion is done 
using Chebyshev polynomials. The equations are time- 
stepped by advancing non-linear and Coriolis terms using 
an explicit second-order Adams-Bashforth scheme while 
the remaining terms are time-advanced using an_implicit 
Crank-Nicolson algor i thm (see iGlatzmaied (jl98j) and 
iChristensen fc Wichtl (|2007f ) for more details). 

3. RESULTS 
3.1. Diagnostic parameters 

In the following discussions we will use several non- 
dimensional diagnostic parameters which are derived 
from the numerical simulations. They describe repre- 
sentative mean values of the flow velocity (or kinetic en- 
ergy) and magnetic field strength (or magnetic energy). 
While in the Boussinesq case these properties have simi- 
lar amplitude throughout the shell, with density stratifi- 
cation their magnitude can vary significantly with radius, 
in particular for the velocity. Hence the question arises 
which property to average in which way. We found that 
averaging energies divided by unit mass gives the best 
results. 

The non-dimensional rms velocity is given by the 
Rossby number 



Ro=JyJ{u-u)dV 



(13) 



where V is volume of the fluid shell and the volume inte- 
gral is evaluated in the fluid shell. The non-dimensional 
kinetic energy per unit mass is then i?o^/2. The Lorentz 
number Lo describes the non-dimensional magnetic fleld 
strength and is deflned as 



Lo = 



I JiB-B)dV 
JpdV 



(14) 



with Ejn — Lo^/2 being the non-dimensional magnetic 
energy per unit mass. The non-dimensional power gener- 
ated per unit mass by thermal buoyancy, scaled by fl^D^ , 
is 

RaE^ J Ur p g s dV 



P 



P. 



JpdV 



(15) 



where u^ is the radial component of u. The Nusselt num- 
ber Nu is the ratio of total transported heat flux to the 
conducted heat flux. The rate of energy dissipated per 
unit mass by ohmic dissipation is 



D. 



The ohmic fraction 



E /A(V X BfdV 



P„. 



JpdV 



fohm — D\IP 



(16) 



(17) 



is the fraction of energy dissipated by joule dissipation 
alone. The characteristic time scale of magnetic energy 
dissipation is 



Err 



(18) 



All diagnostic quantities are time-averaged after a statis- 
tically stationary state has been reached in the simula- 
tions. It would have been more consistent if in Eqs. ()14t - 
[TB)) the magnetic energy, power, and dissipation would 
have been mass-normalized before integration (as is done 
in the case of the kinetic energy) instead of normaliz- 
ing the total energy by the integral of p (which is the 
non-dimensional total mass). However, the integrations 
have been done during the simulation and could not be 
repeated in modified form without repeating the whole 
run. Because the magnitude of these properties seems to 
vary less with radius than that of velocity, the difference 
is probably not critical. 

We separate the dynamos resulting from our numerical 
simulation in two categories: the dipolar category con- 
tains dynamos with a dominant axisymmetric magnetic 
dipole (spherical harmonic degree £ = 1 and harmonic 
order m = 0); the multipolar class contains all other 
kinds of dynamos. To carry out this separation we use 
the dipolarity fdip which is the magnetic energy in the 
axisymmetric dipole normalized by the cumulated mag- 
netic energy in harmonic degrees up to 12, both evalu- 
ated at the outer boundary. Dynamos with fdip > 0.3 
are considered dipolar. It must be noted that our data 
spans the range 10~^ < fdip < 0.9 quite uniformly, hence 
a cutoff of 0.3 remains somewhat arbitrary. We justify 
our choice of cutoff in Sec. 13.41 



3.2. Dynamo database 

Important physical attributes of various numerical se- 
tups used in this study are tabulated in Tab. [TJ The 
database incorporates many important features, such 
as density stratification and variable transport proper- 
ties. The largest density contrast in our simulation is 
Np — 5.5, i.e. Pi/ Po ~ 245. The aspect ratio rj is varied 
to cover dynamos operating in thick shells or nearly full 
spheres as well as dynamos in thinner shells (from r; = 0.1 
to 1] = 0.75). The polytropic index m is also changed in a 
few cases to model different polytropic states. With these 
features our simulations seek to represent, in a simplified 
way, the dynamo action in the fully or partially convec- 
tive and density stratified interiors in different classes of 
objects ranging from giant planets to low-mass stars. 

In total, 272 dynamo cases are simulated. To the 
best of our knowledge this is the largest database of this 
kind. A portion of this database has already been used 



in earlier studies to explore different aspe cts of dynamo 
mecha n ism: A2a, A2b , A5, B la, Bi b fromlGastine et al.l 
(|2012D . iGastineet all (I2013D and iDuarte et al.l (|2013D : 
A3a from lYadav et al.l (|2013[ ). More details, e.g. on the 
choice of model parameters, can be found in these refer- 
ences. The simulations cover a large parameter space: 
1 X 10"'5 < ^ < 1 X 10~3^ 0.2 < Pm < 20, and 
0.3 < Pr < 10. No azimuthal symmetry was imposed 
in the simulations, except for the large aspect ratio run 
of A6 where four fold symmetry in longitude was used. 
All simulations were run for at least one magnetic diffu- 
sion time (D'^/X). 

Our aim is to explore scaling laws in the rotationally 
dominated dynamic regime, i.e. Coriolis forces are domi- 
nant and inertia plays a secondary role, hence we report 
and analyse dynamos with Ixl0^^<£'<lxl0~^ and 
avoid convection with very high supercritical Rayleigh 
numbers. Furthermore, to ensure fully developed con- 
vection in the spherical shell we use only dynamos with 
Nu> 2. 

3.3. Velocity power scaling 

Least-square-optimization of the Rossby number Ro as 
a function of power P and magnetic Prandtl number Pm 
for our database results in 



Ro^2A9- 



P 



0.45 



Pr 



0.13- 



(19) 



This empirical scaling is represented by a solid line in 
Fig. [TJ The dynamo cases marked by a "-I-" in Fig. [T] 
(subset A2b and Bib; see Tab. [Ij are excluded from 
the fit. The excluded cases have exponentially decaying 
electrical conductivity in the outer parts of the spheri- 
cal shell. In these regions Lorentz forces are weak which 
allows for the developm ent of strong z onal flows that en- 
hance the rms velocity (iDuarte et al.l 20131. This is evi- 
dent in Fig. [1] where most of these cases have higher Ro 
and lie above other dynamo cases. The dynamo cases 
from subset A4 with moderate conductivity variation, 
which are marked by a black dot in Fig. [l] do not show 
enhanced zonal flow and their Rossby numbers agree well 
with that of the other dynamos. 

The scatter of data points in Fig. [1] is primarily 
due to differences i n the relative strength of zonal 
flows (jYadav et al.l I2013D . Magnetic field geome- 
try, which can vary substantially from one simulation 
to anot her, aff ects the z onal flow through Maxwell 
stresses (jBrown ing 200^. lYadav et all (|2013D reported 
two different Ro scalings for dipolar and multipolar dy- 
namos (subset A3a; see Tab. [IJ, although the difference 
was relatively minor. They argued that this offset is due 
to more efficient zonal flow quenching by magnetic brak- 
ing in dipolar dynamos as compared to the multipolar 
ones. Such consistent dichotomy is washed out in Fig. [1] 
probably due to our much more diverse da taset. 

On similar lines as lYadav et al.l (|2013[ ). we define a 
convective Rossby number Roconv which is calculated by 
excluding zonal and meridional flow components. A best- 
fit of Roconv gives 



ROconv = 1-66 



pO.- 

pO.08 ' 



(20) 



which is portrayed in Fig. [2l As compared to Fig. [T] the 



Table 1 

Simulation setups 



Subset name 


Core 


BC 


^ 


q(r) 


m 


Np 


Varying inside 


(number of cases) 














the shell 


Anelastic 


















Al(8) 


Insulating 


SF 


0.1 


r 


2 


2, 4 




none 


A2a (44) 


Conducting 


mixed 


0.2 


r 


2 


0.5-5 




none 


A2b (43) 


Conducting 


mixed 


0.2 


r 


2 


1 -5 




A 


A3a (57) 


Insulating 


SF 


0.35 


r 


2 


0.5-5 




none 


A3b (9) 


Insulating 


SF 


0.35 


r-^ 


1.5 


1, 3 




none 


A4 (11) 


Insulating 


SF 


0.35 


r-^ 


2 


3, 5 




U, K, X 


A5 (24) 


Insulating 


SF 


0.6 


r-^ 


2 


0.01 -3 




none 


A6(l) 


Insulating 


SF 


0.75 


r-^ 


2 


1 




none 


Boussinesq 


















Bla (14) 


Conducting 


mixed 


0.2 


r 










none 


Bib (21) 


Conducting 


mixed 


0.2 


r 










A 


B2 (40) 


Insulating 


SF 


0.35 


r 










none 



Notes: BC - mechanical boundary conditions, where SF stands for shear stress free on both boundaries and mixed stands for a rigid inner 
and stress-free outer boundary, aspect ratio rj, gravity profile g(r), polytropic index m, and number of density scale heights inside the shell 
Np for various simulation subsets. Th e la st columns list transport properties which va ry along the shell radius. In groups A2b and Bib 
the variation of A is obtained from Eq. 1111 and in group A4 the variation follows Eq. 1121 
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Figure 1. Rossby number versus a combination of non- 
dimensional power and magnetic Prandtl number. The solid line 
represents Ro = 2.49 P"'*^ Pm""-^-^. Filled (empty) symbols dis- 
play dipolar (multipolar) dynamos. The Ekman number is color 
coded and the symbol shape indicates degree of density stratifi- 
cation. Symbols containing a "-I-" marker have exponentially de- 
caying electrical conductivity in the outer regions (subset A2b and 
Bib in Tab.lTJ and the ones carrying a dot symbol have moderate 
conductivity variations (subset A4 in Tab.[TJ. 

scatter in Fig. [5] is much smaller and also cases with a 
poorly conducting exterior region are now well fitted. On 
the other hand, the Rossby number based on the zonal 
flow component is strongly scattered in our database 
(not shown) and does not exhibit a consistent scaling, 
in contrast to what has been foun d in purely hydro- 



2001 



dynamic spherical shell convecti o n (lAubert et al.l 
CJiristensen' '2002"; 'Aubert] l2005t iShowman et al.l 
Gastine fc Wicht ,2012.). As discussed above, this is due 



2011 



Figure 2. Convective Rossby number versus a combination of 
power and magnetic Prandtl number. The solid line represents 
Roconv = 1.66 PO'*^ P^-OOS, Refer to Fig.[T]for symbol descrip- 
tion. 

to the variable influence of magnetic field on the zonal 
flow via Maxwell stresses. 

A scaling law of the type Ro ex P" with 
0.41 < a < 0.44 has been reported by many 
studies of Boussinesq numerical dynamo simula- 
tions: in a wide control parameter study with rigid 
mechanical boundaries ([Christens en fc AubertI 120061 : 
iTakahashi et al. 20081 ) ; for compositional and thermal 
convection with variable core size, including cases wit h 
nearly fully convective interiors (jAubert et al.l I2009D : 
in si mulations with diff erent thermal boundary condi- 
tions (jChristensenl 12010); and for cases w ith stress-free 
mechanical boundaries (jYadav et alJl2013D . 



The theoretical justifications for the velocity scaling 
Ro ex P", with a slight l y larg er than 0.4, is not en- 
tirely clear. iChristensenI (|20100 discussed several pro- 
posed scaling approaches and showed that they all lead 
to a law of this form, but with differe nt exp onents. The 
mixing length theory ( Bohm-Vitens^ I1958[) applied to 
planetary interiors predicts a = 1/3, but is probably 
not ap plicable to rapidly r otating systems (jChristensenl 
I2O1O0 . lAubert et air()2001f ) considered a triple force bal- 
ance of Coriolis, inertial, and Archimedean forces (the 
so-called CIA balance) in the hydro dynamic fluid vortic- 
ity eq uation and derived a = 2/5. IStarchenko fc Joiiesl 
(|2002l) obtained a = 1/2 based on a balance of Coriolis 
and buoyancy forces. The prediction of the CIA balance 
is the closest to the exponent found by fitting the data 
from dynamo s imulations, but underpredicts it slightly. 
IDavidsoiil ([2013) presented a scaling theory starting from 
the assumption that the magnetic field is independent 
of rotation rate and scales with the cubic root of the 
power, as confirmed by th e numerical results (see next 
section). IDavidsoiil ()2013[) considers a triple force bal- 
ance of Lorentz, Archimedean, and Coriolis forces (the 
so called MAC balance) to then derive a velocity scaling 
of the form Ro ex P^^^. The exponent is very close to 
our empirically obtained exponents a = 0.45 (Eq. ITO|) or 
a — 0.42 (Eq. I^Hl) . None of the scaling theories for ve- 
locity proposed so far accounts for the P„j dependence. 
Inclusion of Pm as a fit parameter signi ficantly reduces 
the scatter of data-points {C hristcnscn fc Auber^ 120061 : 
lYadav et al]|2013HStelzer fc Jackson..2013ir 

3.4. Magnetic field scaling 

In t h e Boussinesq dynamo simul ations. IChristensenI 
(poTol) . ISchrinner et all (|20T1 . and lYadav et al.l (120131 ) 
found that the magnetic field strength obeys very simi- 
lar scaling laws for dipolar and for multipolar dynamos, 
except for a difference in the pre-exponential constant. 
Therefore we treat the two classes of dynamos separately. 
A least-square fit of the Lorentz number corrected by 
fohrn as a function of P and Pm leads to 
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V^ 



= 1.08 P°-^^P, 



km 



for the dipolar dynamos and 
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vfohr, 



0.78 P' 



0.34 



0.13 



0.09 



(21) 



(22) 



for the multipolar ones. These two shifted scalings pro- 
vide a basis for separating dynamos in dipolar and mul- 
tipolar categories. We found that fdip = 0.3 provides a 
reasonably good cutoff; lower or higher cutoff's start to 
mix data points and we loose the demarcation noticeable 
in Fig. 121 

The two expressions (I^TI) and (|22p are very similar ex- 
cept for the larger prefactor in the dipolar scaling. We 
therefore attempt to fit both classes of dynamos with the 
same approximate exponents for P and P,n 



Lo 



vJohr. 



C P 3 P,„ To 



(23) 



where only c is fitted separately for dipolar and for mul- 
tipolar dynamos. This approximate scaling is plotted in 
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Figure 3. Lorentz number corrected for fraction of oliniic dissi- 
pation versus a combination of non-dimensional power and mag- 
netic Prandtl number. The solid line fitting the dipolar dynamos 
(filled symbols) represents Lo/v^/ota^ = 0.91 P^/^ p^l/lO „hile 
the dashed line fitting the multipolar dynamos (open symbols) rep- 
resents Lo/V/oh^ = 0.70 Pl/^ Pm^^^°. Refer to Fig.nfor symbol 
description. 



Fig. 13 lYadav et all (|2013[ ) argued that the higher field 
strengths in dipolar dynamos can be attributed to lower 
fiow velocity in dipolar dynamos (lack of strong zonal 
fiows). Such argument may not work for our database 
since both dipolar and multipolar dynamos are consis- 
tent with a single velocity scaling (see previous section). 
Further study is required to explain this behaviour. 

Two different theoretical approaches have been put 
forward to explain the scaling of the Lorentz num- 
ber with (approximate l y) the cubic root of the power. 
iChristensen fc AubertI ()2006D take a velocity scal- 
ing Ro (X P""^ as given and use the results of 
IChristensen fc Tilgneii (|2004f ) that the magnetic dissi- 
pation time scale Tmag (in rotational units) is approx- 
imately proport ional to Ro~^ . In a different app roach, 
IDavidsonI (|2013') uses the finding bv lKunnen et all (|2010[ ) 
that the small scale vorticity in rotating convection is in- 
dependent of rotation rate. Assuming that this also holds 
for vorticity at the magnetic dissipation length scale Imin , 
for which the local (scale-dependent) magnetic Reynolds 
number is of order one, he infers that also ?mi„/A ex Tmag 
must be independent of the rotation rate. Setting the 
fohrn factor asidc, the magnetic energy is the product 
of power and magnetic dissipation time and is therefore 
also independent of the rotation rate. If magnetic field is 
only a function of power, dimensional arguments dictate 
that it must depend on the cubic root of the power. 

Davidson's scaling theory predicts a dependence of 
the magnetic diffusion time on the Rossby number as 
Tmag oc Ro~^/'^, significantly wea ker than the Ro~^ de- 
pendence originally envisaged bv IChristensen fc Tilgneii 
(|2004) . The scaling of Tmag vs. Ro in numerical dy- 




Figure 4. Magnetic dissipation time versus Rossby number. The 
solid line represents Tmag = 0.61 -Ro""'^^. Refer to Fig. [l] for 
symbol description. 

namo model s had b een r evisited by iChristensed (I2010D . 
lYadav et al] (|2013D . and iStelzer fc Jacksod ()2013f ) who 
aU found exponents somewhat weaker than -1. In Fig. 2] 



we plot 



for our data set against Ro. The scatter is 



larger than in the cases of fitting Ro and Lo, however, 
we find that a dependence 



0.61 i?o" 



(24) 



fits the data best. This scaUng ex ponent is very close to 
the prediction bv iDavidsonI (|2013[ ). 

4. DISCUSSION AND CONCLUSIONS 

In this study we analyze a large number of numerical 
simulations of dynamo in rapidly rotating spherical shells 
to infer empirical scaling laws for mean velocity and mean 
magnetic field. Our database contains anelastic com- 
pressible dynamos with different shell thicknesses, dif- 
ferent density stratifications with density contrast upto 
~ 245, and radially varying transport properties (vis- 
cosity, thermal diffusivity, magnetic diffusivity) . These 
models capture, in a simplified way, conditions that are 
relevant for rapidly rotating gas planets and low-mass 
stars. 

Our study shows that it is possible to relate flow ve- 
locity and magnetic field strength to the available power, 
by a single scaling law in each case, for a wide range of 
different conditions. In particular the same scaling law 
covers incompressible cases and cases with strong den- 
sity stratification. Our finding sancti fies the application 
of the same magnetic field scaling bv iChristensen et al.l 
(|2009| ) to such diverse objects as the Earth's core and 
low-mass stars. 

Our best-fitting power exponent of 0.45 for the Rossby 
number, as a measure of the characteristic velocity, is 
very close to the value of 4/9 obtained by the scaling the- 
ory of iDavidson. (201^) which assumes a MAC balance. 



Our data are better collapsed when the Rossby number 
is calculated only with the convective flow (excluding the 
zonal and meridional flow component) and the optimal 
exponent of 0.42 is intermediate between Davidson's pre- 
diction and the valu e of 2/5 inferred from a CIA balance 
(jAubert et al.l 120011 ). Because our simulations may not 
have exactly reached an asymptotic regime, the precise 
value of the exponent and its theoretical justification re- 
main somewhat uncertain. 

Our data for the Lorentz number which measures mag- 
netic field strength are fully compatible with a cubic root 
dependence on the non-dimensional power (after correct- 
ing for that part of the energy flux that does not con- 
tribute to magnetic field generation). This makes the 
scaling independent of the rotat ion rate, which i s a cen- 
tral point in the scaling theory ofl DavidsonI (|2013| ) . While 
both dipolar and multipolar dynamos follow basically the 
same scaling law, we find that in t he dipolar ca s e the field 
strength is larger by « 30%. In lYadav et al.l (|2013D we 
attributed this difference to a slightly weaker flow veloc- 
ity in dipolar dynamos, but in our present study we do 
not find a significant difference in the Rossby number 
between dipolar and multipolar cases. However, we note 
that the ohmic dissipation time is systematically shorter 
for multipolar dynamos (open symbols in Fig. U]) than it 
is for dipolar dynamos (filled symbols), i.e., the gener- 
ated magnetic energy is more rapidly dissipated in the 
multipolar case, leading to a lower mean field strength. 
The reason for the longer dissipation time of dipolar dy- 
namos could be that more magnetic energy is stored in 
large spatial scales where it is less vulnerable to dissipa- 
tive destruction. 

There are few outstanding issues which need fur- 
ther explorati on. As ha s recently been emphasized by 
IStelzer & Jackson! (|2013( ). some dependence on the mag- 
netic Prandtl number is required for an optimal fit of 
the dynamo simulation data. Although the influence 
of Pm in the scaling laws is secondary, its value in 
natural dynamos is several orders of magnitude lower 
than it is in the simulations. Ignoring or including 
Prn in the scaling laws makes a difference of nearly 
an or der of magnitude in the pred iction for planets or 
stars. IChristensen fc Tilgneii ()2004[ ) argued that the Pm- 
dependence, although needed for fitting the model data 
at Pm, ~ 1, would disapp ear in the astrophysica l ly rel - 
evant range Pm ^ 1 and IChristensen fc Aube rt' (2006|) 
found good agreement between prediction and observa- 
tion in case of the Earth's dynamo only when the depen- 
dence on the magnetic Prandtl number is ignored. How- 
ever, the assumption that the Pm-dependence vanishes 
at small values has not yet been confirmed by theory or 
by numerical simulations, which are not available for P,„ 
much smaller than one. 

Another issue is the independence of magnetic field 
scaling from rotation rate that is implied by our scaling 
law. For low-mass stars a saturated (rotation-rate inde- 
pendent) magnetic surface flux is observed for very rapid 
rotators with low Rossby number, but for larger values 
of the Rossby number it decreases r oughly proportionall y 
with decr e asing rotation rate (e.g. iReiners et all 120091) . 
iSchrinneil ()2013l ) argues that the rotational dependence 
at large Rossby number may enter through its influence 
on the ohmic dissipation factor fohrm challenging the of- 
ten made assumption that fohm ~ 1 in all kinds of nat- 



ural dynamos. 

One problem with applying our scaling laws to stars 
and planets is that they predict average values inside 
the dynamo, whereas observations usually relate to the 
strength of the magnetic field or part of the field at the 
surface of the object. The assumption that the two values 
can be related by a fixed factor is probably too simplis- 
tic. It is conceivable that a dependence of this factor on 
the rotation rate could contribute to the observed depen- 
dence of the surface magnetic flux on the Rossby number 
in M-stars with moderate to low rotation rates. 

Our work extends earlier stud- 

ies (iChristensen fc Aube^ 120061: pTakahashi et all 120081: 



Aubert et all 120091: iChristensenl EolOl: iSchrinner et all 
20121 : lYadav et alIl2013D . It shows that previously found 



scaling laws for velocity and magnetic field strength in 
rotating dynamos also apply when strong radial varia- 
tions in density or transport properties are present, as 
expected in the conducting cores of gas planets, brown 
dwarfs and low-mass stars. This strongly supports the 
notion that the power-based scaling laws are rather 
universal. For the magnetic field scaling there is a decent 
agreement with observation for sev eral solar system 
planets and certain cla sses of stars (jChristensen et aTl 
120091 : [Christenseiill2010D . However, the velocity scaling 
has not yet been tested against observation (except for 
Earth's core). This together with the open questions on 
the role of the magnetic Prandtl number, the cause for 
weaker surface magnetic fields at less rapidly rotating 
stars, and the relation between internal and surface 
magnetic fields keep this research area vibrant and 
provide an exciting avenue for future explorations. 
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(HLRN), and the Western Canada Research Grid 
(WestGrid) and Reseau quebecois de calcul haute 
performance (RQCHP). All the figures were generated 
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